Insights from ARPES for an undoped, four-layered, two-gap high-T c superconductor 
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An undoped cuprate with apical fluorine and inner (i) and outer (o) CuCh-layers is a 60 K super- 
conductor whose Fermi surface (FS) has large n and p-doped sheets with the SC gap on the n-sheet 
twice that on the p-sheet Q]. The Fermi surface is not reproduced by the LDA, but the screening 
must be substantially reduced due to electronic correlations, and oxygen in the o-layers must be al- 
lowed to dimple outwards. This charges the i-layers by 0.01 1 e| , causes an 0.4 eV Madelung-potential 
difference between the i and o-layers, quenches the i-o hopping, and localizes the n-sheets onto the 
i-layers, thus protecting their d-wave pairs from being broken by scattering on impurities in the BaF 
layers. The correlation-reduced screening strengthens the coupling to z-axis phonons. 

PACS numbers: 71. 18.+y, 74.25. Jb,74.72.Jt 



Despite 20 years' intensive research, high-temperature 
superconductivity in the cuprates has not been under- 
stood. Much insight has been gained through discoveries 
of new cuprates and refinements of experimental tech- 
niques. An example is the discovery |2| of an undoped 
4-layered cuprate with apical fluorine, Ba2Ca3Cu40sF2, 
which is not a Mott insulator, but a 60 K HTSC; the ex- 
perimental findings about its Fermi surface and SC gap 
are reported in the preceding Letter 0. 

All HTSCs consist of Cu0 2 -layers separated by api- 
cal blocks into which dopants may be inserted. Un- 
doped cuprates (Cu 3d 9 ) are antiferromagnetic (AF) 
Mott insulators. Upon hole-doping beyond p~0.05 
per Cu02, they become HTSCs with pairing sym- 
metry d x 2_ y 2, and the critical temperature reaches a 
maximum, T cmax , when p~0.2. For optimally and 
over-doped materials the Fermi surfaces (FS) measured 
by ARPES enclose large hole-areas of size 1 + p (in 
units of half the Brillouin-zone area), are centered at 
(k x , k y ) = (tt, 7t) , and agree surprisingly well with de- 
tailed predictions from LDA band-structure calculations 
3J. Even the existence of a small sub-band splitting 
in bi-layered cuprates 0, , which is twice the integral 
for inter-layer hopping, t± (k„) ~ t± (cos k x — cos k y ) 2 /4, 
and the existence of a k z -dispersion proportional to 
t± (k M ) cos \k x cos \k y cos ^k z in body-centered tetrago- 
nal materials 0, 0], have recently been confirmed by 
ARPES [3J, |S| and angular magnetoresistance oscilla- 
tions (AMRO) 0. Hence, over- and optimally doped 
cuprates appear to have well-defined quasi-particle bands 
with mass-renormalizations 2-4, not only in the direc- 
tion of the layers, but also in the perpendicular direction 
I n the superconducting state, the magnitude 
of the gap is approximately 5.bksT c . For under-doped 
materials, a pseudogap opens up when T < T*. The re- 
maining Fermi arcs |3|, Il2| . when mirrored around the 
AF zone boundary, enclose small hole-areas of size p, 
are centered at (f , f ) , and are roughly reproduced by 
LDA+f/ calculations with Udd^QeV and AF long-range 



order, or similar Hartree-Fock (HF) solutions for 3-band 
models 0, 0|. Upon doping with electrons beyond a 
certain level, a metallic AF (AFM) phase develops and, 
subsequently, a superconducting phase with the electrons 
occupying small, (tt. 0)-centered FS areas of size n. For 
higher n-doping |3|, Il4j . antiferromagnetism disappears 
and the FS again becomes (tt, 7r)-centered and large, of 
size 1— n. Optimal n-doping occurs for n~0.2, and known 
values of T c max are much smaller (< 40 K) than for opti- 
mal p-doping (< 140 K) , and so are the gaps. So far it has 
not been possible to p- and n-dope the same compound. 

In Fig.^ we show the FS of undoped Ba2Ca3Cu408F 2 . 
ARPES resolves the two sheets shown in grey Q]. Both 
are (tt, 7r)-centered and large, their areas being l±p with 
p~0.2. We shall refer to the larger one as the hole-doped 
or p-sheet and to the smaller one as the electron-doped 
or n-sheet. Ba 2 Ca3Cu40sF2 thus appears to be the first 
se//-doped HTSC having n- and p-doping in one and 
the same crystal. All other known undoped cuprates, 
single- as well as multi-layered, are insulating, and when 
hole-doped, the observed FS-area splittings |8| are pro- 
portional to t± (k M ) and almost an order of magnitude 
smaller than the ±0.2 found in Ba 2 Ca3Cu40sF2. The 
question we shall try to answer in this letter is: How and 
why does this compound manage to dope itself? 

There are no FS measurements for multi-layered 
cuprates with mequivalent inner (i) and outer (o) layers, 
with the exception of Ba2Ca3Cu40sF2, but site-selective 
NMR studies indicate that i-layers have less holes 
than o-layers and that superconductivity develops dif- 
ferently with temperature in the different layers. This 
behaviour is interpreted in terms of electrons localized 
onto individual layers, each of which follows the generic 
phase diagram as a function of doping, and which are 
weakly coupled by proximity effects [Tg. Most spectac- 
ular are the properties of 5-layered HgBa2Ca4Cu50i 2+( 5 
[T7^ . Here, the o-layers form an optimally doped HTSC 
with T c max =108K and the three i-layers are AF metals 
(AFM) with moments of about 0.3/is. By decreasing the 
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FIG. 1: FS for BaaCaaCruOgFa. Grey: ARPES [lj. The 
width reflects the experimental error bar. The smaller (larger) 
sheet is the n (p)-sheet. Open dots: Calculations. Top right: 
LDA. Remaining three quarters: (LDA+(/). Bottom right: 
Flat CuO-2-layers and 2% F-0 exchange (g=0.02) treated in 
the VCA. Indistinguishable from this FS are the ones calcu- 
lated with no F-O exchange, but with O in the o- or i-layers 
dimpled outwards by 0.07 A. See Fig.|H] The n (p)-sheet re- 
sides on the i (o)-layers. ef is adjusted such the FS fits along 
the nodal line whereby n=0.19 andp=0.12. Bottom left: Dim- 
pling as before, but with the vertical Ba and F positions ad- 
justed such as to give the best agreement with the ARPES 
FS: «Ba — £Cu is decreased from 2.0 to 1.7A and zf — zcu Q is 
decreased from 2.4 to 2.lA. ef is adjusted as before, whereby 
n=0.18 and p=0.21. Top left: As before but with O in the 
o- (or i-) layers dimpled inwards by 0.10 A. Now, the n(p)- 
sheet resides on the o(i)-layers. ZBa — zcu = 1.6A and 
zf — zcu a = 2.0A. ef is adjusted as before, whereby n=0.15 
and p=0.18. Note that inter-sheet splittings in k-space are 
enhanced near the saddlepoints in the band structure. 

hole doping, T c decreases to 70 K and the o-layers form a 
uniformly mixed HTSC-AFM phase with O.l^ig moments 
while the i-layers form a 0.7u£-AFM. The layers are flat 
and the estimated hole-counts for the o- and i-layers dif- 
fer by as much as 0.2 for the over-doped material. 

Seen on the background of these truly exotic proper- 
ties, the discovery |l| that Ba2Ca3Cu40sF2 is the first 
HTSC with different SC gaps on the two sheets, may 
seem less spectacular. However, it is the gap on the n- 
sheet which is the largest, actually twice the gap on the 
p-sheet, and this is opposite to what one would expect 
from the generic phase diagram. Moreover, two-gap su- 
perconductivity has only been observed unambiguously 
in a few materials of which MgB2 is the mo st sp ectacular. 
There, the mechanism is conventional 

Whereas the FS sheets observed in all other multi- 
layered HTSCs are hardly split along the nodal direction, 



Ba2Ca3Cu40sF2 exhibits a large nodal splitting, Afcmoi- 
This is the third anomaly of the ARPES data and the 
key to the self-doping and the two gaps. For compari- 
son with the ARPES FS, we show in the top right-hand 
quarter of Fig^the LDA FS which is seen to have four 
FS sheets, split almost exclusively by inter-layer hop- 
ping [2(| • The four sub-bands are split by approximately 
±i (y/5 ± l) t±_ (k„) with t_L~0.2 eV. This is the first case 
known to us where there is a substantial discrepancy be- 
tween an experimental, large FS and the LDA. The lat- 
ter fails to reproduce the Afcrno] de/dku\Q\ = 0.4 eV large 
Madelung-potential difference (on electron scale) between 
the layers on which respectively the p- and the n-nodal 
quasi-particles are located. We have taken Qe/dfcnio] as 
the bare (LDA) nodal velocit y, w hich is the consistent 
choice in the present context [21|. In the LDA (GGA), 
the electron potential averaged over an i-layer is 0.024 eV 
(0.034 eV) less than that averaged over an o-layer, that is 
17 (12) times too small! The discrepancy remains after 
we optimize the structure using the LDA. Thereby, most 
notably, the oxygens in the o-layers dimple outwards by 
~3°, about one third the dimple in YBa 2 Cu307 |2^ . 

Since the 0.4 eV potential exceeds \ (v5 + l) t±, it ef- 
fectively blocks the hopping between the i- and the o- 
layers. Therefore, the n-sheets are associated with the 
i-layers and the p-sheets with the o-layers, or the other 
way around. The potential difference thus blocks inter- 
band scattering and protects the «-bands from scattering 
on the ubiquitous impurities in the apical BaF-blocks. 
This may be the reason why two SC gaps survive and 
-in case the p-sheets reside on the o-layers- also the rea- 
son why the gap on the p-sheets is the smaller one: it is 
suppressed by scattering on impurities in the BaF-block. 

Let us set this 0.4 eV potential difference in perspec- 
tive: With the in-plane lattice constant 3.86 A and the 
distance between Ca and Cu02 layers 1.54A, Poisson's 
equation says that the layer-averaged potential is a saw- 
tooth with kinks of ^20eV<5 ra on the layers. Q n is the 
charge per 2D cell. If Ba2Ca3Cu40sF2 were purely ionic, 
the potential would therefore be on the Ca 2+ layers and 
20 eV on the (CUO2) 2- layers, the same on the i- and 
o-layers because the CaCu02-wirf is neutral. Next, we 
transfer from the BaF layers qi positive charges to each of 
the «-layers and q positive charges to each of the o-layers. 
This sets up a difference in Madelung potential between 
o- and i-layers of 40g^ eV, which depends only on qi, not 
on q a .The ARPES value of 0.4 eV (neglecting exchange- 
correlation contributions) therefore tells us that the i- 
layers have a net charge of either 0.01 |e|, in which case 
the n-sheets reside on the i-layers, or of — 0.01 |e|, in 
which case the p-sheets are on the i-layers. In the LDA 
(GGA), the charging is 17 (12) times smaller. 

Crystallographically it is possible that q of the two 
layer-oxygens per cell are exchanged with fluorine, so 
that the formula is Ba2Ca3(Cu02- g F g )404 g F2_4 g . In 
fact, LDA calculations with g=0.25, employing not only 
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FIG. 2: Total energy as a function of the Raman-active dim- 
pling displacements calculated with LDA+LX See Fig. [3] 
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FIG. 3: Layer-averaged potentials calculated with (LDA+C/) 
for dimplings of ±0.1 A of the inner (left) and outer (right) 
layer-oxygens. This potential is the calculated Cu and O core 
levels, lined up on the flat layer and normalized to zero at the 
Ca mirror plane, m is the mirror plane. 



the virtual-crystal approximation (VCA), but also super- 
cells, do reproduce the ARPES FS. However, the posi- 
tion of the apical F 2p z level is 4 eV below that of apical 
O 2p z , and such impurities in the Cu02-layers are very ef- 
ficient in breaking d a ,2_ J( 2-wave pairs: For an allowed T c - 
suppression of 10%, we find that q cannot exceed 0.001, 
which is even smaller than the charging of 0.01 needed 
to create the observed crystal-field splitting in the ab- 
sence of screening [23J. This rules out F-0 exchange as 
an explanation for the anomalous ARPES FS. 

QMC calculations for a ii„-band Hubbard model for 
Ceo have revealed that the RPA describes metallic sceen- 
ing well for £//W<2, but that once E//W~3, electronic 
correlations reduce the screening by an order of mag- 
nitude [iij]. This seems to be the problem with the 
LDA when applied to Ba2Ca3Cu40sF2, which is even 
nearly 2D. One might think of reducing all inter-layer 
hoppings as was recently done using the Gutzwiller ap- 
proximation |25| , but this would not directly change the 
metallic inter-layer screening, 1 + 2NM ~ 25, where 
2A/~0.6 electrons/ (eV-Cu02) is the intra-layer density of 
states at the Fermi level and M~40 eV-CuO 2 /electron is 
the inter-layer Madelung constant. Instead, we perform 
selfconsistent LDA+C/ calculations for AF (ir, 7r)-order 
and subsequently calculate band-structures and FSs for 
the charge potential, neglecting the spin potential. This 
procedure we denote (LDA+J7). The result for the pure 
compound (q=0) is very similar to that obtained by the 
LDA and shown in the top right-hand quarter of Fig.^ 
However, with an F-0 exchange corresponding to q=0.02, 
we do obtain the proper splitting along the nodal line. 
This is shown in the bottom right-hand quarter of Fig.^ 
The (LDA+Z7) procedure thus reduces the screening of 
this perturbation to about 2, but the pair-breaking ar- 
gument why F-0 exchange is an unlikely explanation for 
the anomalous ARPES FS remains valid. So from now 
on, we take q=0. 

Next, we include the calculated outwards dimple of the 
outer 02-layer (see Fig.|3J), which is about the same in the 
LDA and the LDA+C/. That this may induce a potential 
shift between the i- and o-layers can be seen as follows: 



For the calculated 0.07 A ~ A dimple, the bare value 
of the Cu 2 +-(0 2 -) 2 dipole in the o-layer is 3x2 ,° cV =3eV, 
and after averaging over the Cu and O2 sub-layers this 
yields a 2 eV bare Madelung shift of the o-layer with re- 
spect to the i-layer. Fig-EJ shows the (LDA+C/) layer- 
averaged potential for all the Raman-active (even) dim- 
ples of the i and o-layers. We see that the intra-layer 
screening is substantial, that the inter-layer screening is 
about 4, and that for the calculated equilibrium struc- 
ture, the (LDA+C/) creates the proper 0.4 eV potential! 
The FS is very similar to the one calculated by adjust- 
ing the F-0 exchange to reproduce the nodal splitting 
of the ARPES FS shown in the bottom right quarter of 
Fig-E Essentially the same FS is obtained by dimpling 
the i-layers outwards. In all cases, the n-sheets are on 
the i-layers and split by the hopping between them. The 
p-sheets are on the o-layers and almost degenerate; they 
agree less well with ARPES near the anti-nodal direction 
where the ARPES peaks are less well defined. Moreover, 
the ARPES data has been forced to yield the area, 1 ± p, 
proper for zero doping, which for the real material may 
not hold exactly. We have chosen to shift the Fermi level 
to fit ARPES along the nodal direction. In the bottom 
left-hand quarter of Fig.^ we demonstrate that perfect 
agreement with ARPES can be obtained by adjustment 
of the vertical Ba and F positions. 

If -in disagreement with the total-energy result- we 
dimple the o-layers or the i-layers inwards by 0.1 A, the 
FS shown in the top left-hand quarter of Fig.^ is ob- 
tained. Now the n-sheets are on the o-layers and unsplit, 
whereas the p-sheets are on the i-layers and split by the 
hopping between them. 

There has recently been interest in possible conse- 
quences of poor z-axis screening for the electron-phonon 
interaction at small momentum transfers [26L [23. We 
note that the coupling of the FS to an a\ g Raman-active 
dimpling mode via the screened Madelung-potential dif- 
ference (once it exceeds the inter-layer hopping) is pro- 
portional to the CU-O2 distance, rather than to this dis- 
tance squared, as is the case for the hopping-induced cou- 
pling considered earlier [2^]. From Fig.|21we realize that 
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the deformation potential is sizeable: £>=5eV/A. For a 
rough estimate of the associated electron-phonon cou- 
pling, let us like in Ref. 0] use the Hopfield expression, 
A = ND 2 1 (Mw 2 ) , to compare with MgB 2 for which 
A«1.0 : In Ba 2 Ca3Cu408F 2 , the density of states appro- 
priate for ei x 2_ y 2 -pairing is twice the value in MgB 2 , D 
is 2.5 times smaller, and Mio 2 is 1.3 times smaller for a 
Raman-active oxygen dimpling mode than for the optical 
bond-streching mode in MgB 2 . Finally, the strong screen- 
ing for large momentum transfers in Ba 2 Ca3Cu40sF 2 ef- 
fectively limits the integrand to the central half of the 
2D Brillouin zone. As a result, A~|, which suffices to 
enhance some other mechanism for HTSC. 

We have shown how undoped Ba 2 Ca 3 Cu40gF 2 man- 
ages to self-dope: By a dimpling distortion which is 
screened out poorly due to strong electronic correlations. 
This sets up an 0.4 eV difference in Madelung poten- 
tial between the outer and inner layers, and that causes 
the FS to split along the nodal direction. The large 
Madelung-potential difference prevents interband scat- 
tering, and that seems to be the reason why two SC gaps 
survive. Since the dimpling is outwards according to our 
total-energy LAPW calculations, both in the LDA and 
in the LDA+L 7 , the FS sheet with the small gap, the 
p-sheet, is located on the outer layers, next to the BaF- 
blocks where impurities are ubiquitous. This suggests 
that the p-gap is small because it is suppressed by im- 
purity scattering, while the n-sheet is protected by the 
Madelung potential. The phase-diagram model, on the 
other hand, predicts the p-gap to be the largest. A fur- 
ther reason for the different gap sizes could be that, like 
in MgB 2 [l^.ll9|. (one of) the interaction(s) driving the 
superconductivity is stronger for the n- than for the p- 
sheet. It could be of relevance that the bonding n-sheet 
has a larger effective intra-layer second-nearest neighbor 
Cu-Cu hopping integral, t' + ~ 0.39£, than the p- 
sheet, for which the effective value is merely t' ~ 0.33i, 
and therefore nests better along the [10]-directions. This 
would be consistent with the observation that for hole- 
doped multi-layer cuprates T c max correlates positively 
with the largest sheet- value of t'/t 0. 

Why these anomalies have been observed only in 
Ba 2 Ca3Cu40sF 2 is undoubtedly due to the apical ions 
being F rather than O. Presumably, when going from api- 
cal O to the more ionic F, the covalency between Ba and 
the apical ion is reduced, whereby the covalency between 
Ba and O in the o-layer is strenghtened and causes it to 
dimple towards Ba. This is supported by the experimen- 
tal finding that the outer layers in HgBa 2 Ca4 Cur O i 2+a 
are dimpled inwards [29| or maybe flat [bj, |30j . 

It will be interesting to follow the development of the 
Ba 2 Ca3Cu408(Fi_ 2p 2 p) 2 FS and its SC gaps as a func- 
tion of hole-doping p. If our assignment of the p-sheets 
with the o-layers is correct, and as long as screening re- 
mains poor, O doping should suppress the gap on the 
p-sheet. Currently it is known that T cmax =105K is 



reached at p=0.20 |2J, but no ARPES data are available. 

We are grateful to O. Gunnarsson and J. Zaanen for 
useful discussions, to A. Yamasaki for checking calcula- 
tions. 
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